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Abstract 

A celebrated problem in numerical analysis is to consider Brownian motion orig- 
inating at the centre of a 10 x 1 rectangle, and to evaluate the ratio of probabilities 
of a Brownian path hitting the short ends of the rectangle before hitting one of the 
long sides. For Brownian motion this probability can be calculated exactly [2]. Here 
we consider instead the more difficult problem of a self-avoiding walk in the scaling 
limit, and pose the same question. Assuming that the scaling limit of SAW is confor- 
mally invariant, we evaluate, asymptotically, the same ratio of probabilities. For the 
SAW case we find the probability ratio is approximately 200 times greater than for 
Brownian motion. 



1 Introduction 

In the January 2002 issue of SIAM News L N Trefethen [18J presented ten problems used 
in teaching numerical analysis at Oxford University. The answer to each problem was a 
real number, and readers were set the challenge of computing 10 significant digits of the 
answer to each problem, with the additional incentive of $100 for the best entrant. In 
the event, some 94 teams from 25 nations submitted solutions, and twenty of these teams 
achieved a perfect score. A donor, W J Browning, kindly provided additional funds so 
that all 20 teams achieving a perfect score could be awarded the $100 prize! 

As is so often the case, the mathematics revealed by the different approaches to the 
solutions was at least as interesting as the original problems, and indeed spawned a book 
The Siam 100-Digit Challenge: A study in high-accuracy Numerical Computing, [2j, which 
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is usefully and engagingly described in J Borwein's review [3]. Indeed, the authors of the 
book responded by providing not 10, but 10,000-digit solutions to all the problems but 
one, for which "only" several hundred digits are available. Borwein's personal favourite, 
and the catalyst for the calculation reported here, is problem number 10: A particle at 
the centre of a 10 x 1 rectangle undergoes Brownian motion (i.e., 2-D random walk with 
infinitesimal step lengths) till it hits the boundary. What is the probability that it hits at 
one of the ends rather than at one of the sides? 

Remarkably, it turns out [2] that a closed form solution to this question can be found, 
via the theory of elliptic functions and invoking some results of Ramanujan on singular 
moduli. The result is 

p e = - arcsin ((3 - 2 v / 2) 2 (2 + V5) 2 (VlO - 3) 2 (5 1/4 - V2) 4 ) = 0.000000383758797925 .... 

In discussion with Nathan Clisby, the question as to the corresponding problem where 
Brownian motion is replaced by the scaling limit of a self-avoiding walk (SAW) arose. The 
rest of this paper addresses our answer to this question. In fact, we prefer to address the 
equivalent question, which is the ratio of the probability that the walk hits the end to the 
probability that it hits a side. This is just 

Pe 
l-Pe 

which has the numerical value 0.00000038375894519599411176841999126970034234598936 . . . 
for the random walk case. 

2 Self-avoiding walks 

A self-avoiding walk (SAW) of length n on a periodic graph or lattice C is a sequence of 
distinct vertices wq, w\, . . . , w n in C such that each vertex is a nearest neighbour of its 
predecessor. In Figure [T] a short SAW on the square lattice is shown, 



Figure 1: A self-avoiding walk on the square lattice. 

while in Figure [2] a rather long walk of 2 25 steps is shown (generated by a Monte Carlo 
algorithm [4J). 

One of the most important questions one asks about SAWs is: how many SAWs are 
there of length n, (typically defined up to translations) denoted c n ? Frequently one rather 
considers the associated generating function 

C(x) = Y,CnX n . 
n>0 
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Figure 2: A typical two-dimensional SAW of 2 steps on the square lattice - courtesy of 
Nathan Clisby. 



For SAWs on a lattice it is straightforward to show that c n grows exponentially with n 
[TT] . so that c n oc fi n , where fi = l/x c depends on the lattice, and determines the radius 
of convergence of the generating function C(x). Another important question, and central 
to our investigation, is the existence and nature of the scaling limit of SAWs, described in 
the next section. 

3 The scaling limit 

An intuitive grasp of this concept can be gained by looking at the two figures above. In 
the first, the effect of the lattice is clear. In the second, there is no obvious lattice, and 
indeed no way to tell that this is not a curve in the continuum. We formalise this notion 
as follows: Consider a smooth (enough) closed domain f2, on an underlying square grid, 
with grid spacing 5 as shown in Figure [3j Denote by fts that portion of the grid contained 
in 0. Take two distinct points on the boundary of Q labelled a and b. Now take the nearest 
lattice vertex to a, and label it as, and similarly let b$ be the label of the nearest lattice 
vertex to point b. Consider the set of SAWs u)(fl$) on the finite domain £ls from as to bg. 
Recall that 5 > sets the scale of the grid. Now let \u>\ be the length of a walk ujs £ oj(Qs), 
and weight the walk by x'^' . The reason for this is that the walks are of different lengths, 
making the uniform measure not particularly natural. (There is also a normalising factor, 
which for simplicity we ignore). 

As we let 5 — > we expect the behaviour of the walk to depend on the value of x. For 
x < x c it is possible to prove that us goes to a straight line as S — > 0. (Strictly speaking 
it converges in distribution to a straight line, with fluctuations 0(\/<5).) For x > x c it is 
expected that us becomes (again, in distribution) space-filling as 5 — > 0. But at x = x c 
it is conjectured that ujs becomes (in distribution) a random continuous curve, and is 
conformally invariant. This describes the scaling limit. If this conjecture is correct, a 
second, pivotal, conjecture by Lawler, Schramm and Werner [9] is that this random curve 
converges to SLE g / 3 (described below) from a to 6 in the domain CI. 

3.1 Schramm Lowner Evolution 



For an approachable discussion of SLE K , see Chapter 15 of [6j. Here we give a very 
minimal outline. Let HI denote the upper half-plane. Consider a path 7 starting at the 
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Figure 3: Discretisation of domain f2. 



boundary and finishing at an internal vertex. Then ET\7 is the complement of this path, 
and is a slit upper half-plane. It follows from the Riemann Mapping Theorem that it 
can be conformally mapped to the upper half-plane. Lowner |10j discovered that by 
specifying the map so that it approaches the identity at infinity, the conformal map so 
described (actually a family of maps, appropriate to each point on the curve) satisfies a 
simple differential equation, called the Lowner equation. The mapping can alternatively 
be defined by a real function. This observation led Schramm to apply the Lowner equation 
to a conformally invariant measure for planar curves. That is to say, the Lowner equation 
generates a set of conformal maps, driven by a continuous real-valued function. Schramm's 
profound insight was to use Brownian motion Bt as the driving function^] So let Bt, t > 
be standard Brownian motion on M and let k be a real parameter. Then SLE K is the 
family of conformal maps gt : t > defined by the Lowner equation 

8 ( \ 2 < \ 

dt g t (z) - ^jKB t 

For k < 4 the domain of gt, which is the set of initial conditions z for which the above 
differential equation still has a solution at time t, equals the half-plane minus a simple 
curve. The curve starts at the origin and can be shown to go to oo. This is called chordal 
SLE K as it describes paths growing from the boundary point and ending at the boundary 
point oo. Larger values of k lead to more complicated behaviour. 

For a general simply connected domain D, SLE K between two boundary points is just 
defined by taking the image of SLE K in the half-plane using a conformal map that takes the 
half-plane to D, and and oo to the two prescribed boundary points. It is conjectured, and 
widely believed, that the scaling limit of the SAW between two fixed points is conformally 
invariant. Lawler, Schramm and Werner [9] showed that this conformal invariance and 
the so called restriction property for the SAW imply that the scaling limit of the SAW is 
described by SLE 8 / 3 . Note that this only describes the SAW in D between two boundary 

1 lt is the only process compatible with both conformal invariance and the so-called domain Markov 
property. 
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points. If we want to describe the SAW which starts at a fixed point on the boundary but 
can end anywhere on the boundary, then we need to know the distribution of this random 
endpoint. We discuss this distribution in the following section. 

Hopefully this rather vague description will convey the flavour of this exciting and 
powerful development in studying not just two-dimensional SAWs, but a variety of other 
processes, such as percolation, the random cluster model, and the Ising model. We refer 
the reader to pQ for greater detail of both SLE and these applications. 

3.2 Conformal mapping 

Let D be a bounded, simply connected domain in the complex plane containing 0. We 
are interested in paths in D starting at and ending on the boundary of the domain. 
Initially we will consider random walks, later self-avoiding walks. As described above in 
the description of the scaling limit, we can discretize the space with a lattice of lattice 
spacing 5. In both the random walk case and the SAW case we are interested in the scaling 
limit as 5 — > 0. For random walks the scaling limit is Brownian motion, stopping when 
it hits the boundary of D. The distribution of the end-point is harmonic measure. If the 
domain boundary is piecewise smooth, then harmonic measure is absolutely continuous 
with respect to arc length along the boundary [12]. Let Hd{z) denote the density with 
respect to arc length (often called the Poisson kernel) . If / is a conformal map on D that 
fixes the origin, and such that the boundary of f(D) is also piecewise smooth, then the 
conformal invariance of Brownian motion implies that the density for harmonic measure 
on the boundary of f(D) is related to the boundary of D by 



In the case of the SAW, Lawler Schramm and Werner [9] predicted that the corresponding 
density of the probability measure p(z) transforms under conformal maps as 



where 6 = 5/8 and the constant c is required to ensure that pd(z) is a probability density. 
If one starts the random walk or the SAW at the center of a disc, then the hitting density 
on the circle will be uniform. So the above equations determine the hitting density for 
any simply connected domain. Simulations in [5], [Z] and [8] provide strong support for 
the conjectured behaviour. (Eq. Q is correct for domains whose boundary consists only 
of vertical and horizontal line segments. For general domains there is a lattice effect that 
persists in the scaling limit that produces a factor that depends on the angle of the tangent 
to the boundary that must be included [7j.) 

The solution of the original problem for random walks by conformal maps is described 
in P|, where a conformal map from the unit disc to an a x c rectangle is given by the 
Schwarz-Christoffel formula. Our approach also uses conformal maps, but we instead use 
a map between the upper half-plane and a rectangle, where the mapping is again given by 
a Schwarz-Christoffel transformation. For a > 1, let 



h D (z) = \f'(z)\h f{D) (f(z)). 



(1) 



p D (z)=c\f(z)\ b p f{D) (f(z)) 



(2) 
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f(z) is a Schwarz-Christoffel transformation that maps the upper half plane to a rectangle. 
The rectangle has one edge along the real axis and is a midpoint of this side. So the 
corners can be written as ±a/2 and ic±a/2 where a, c > are the length of the horizontal 
and vertical edges, respectively. We have 

/(l) = a/2, /(-l) = -a/2, /(a) = a/2 + ic, /(-a) = -a/2 + ic, /(0) = 0. 

So 

a 



dx f a dx 



We note that 



-i vT — x 2 V« 2 — £ 2 A V x 2 — \\J a 2 — x 



-K | - 1 . r - — K 1 — I . n > I . 



a \a 




Here K(x) is the complete elliptic integral of the first kind. By dilation invariance we only 
need concern ourselves with the aspect ratio a/c. So given an aspect ratio r, we have to 
find q such that 

J\{l-x 2 r l l 2 (a 2 -x 2 )- l l 2 dx 2K(i) 



~ f"(x 2 - 1)-V2( a 2 _ x 2)-i/2 dx K ^ v^ET^j ' 

Setting k = 1/a, then c can be written K(/c') = K'(A;). Therefore 

c _ 1 _ K(l/a') _ K'(l/a) 
a~ r ~ 2K(l/a) ~ 2K(l/a) " 

Now the nome g is defineej^Jby g = exp(— 7rK'/K), so 

_ 2K(l/a) 
r ~ K'(1/q) ' 

and 

e- 27r / r = g. 

Then we have 
so 

- ^ 3 (e- 2 ^) 



" 2 (e- 2 ^)' 

where 0j(g) = 0j(0,q) is the Jacobi theta function. Evaluating this in one's favourite 
algebraic package gives the required value of a for any r > 1 instantly. 

Alternatively, we note that for an aspect ratio of 10, as posed in the original problem, 
a is very close to 1. We can achieve very high accuracy by expanding the ratio of the 
above integrals around a = 1 from the Taylor expansion of elliptic integrals, and find 

r = i (41og(2V2) - 21og(« -!) + («-!)- |( Q - l) 2 + A (a - l) 3 + 0(a - l) 4 ) . 



2 Amramowitz and Stegun, p591 formula 17.3.17 
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Solving this numerically for r = 10 gives a = 1.00000120561454706472212 .... To leading 
order one obtains 

a « 1 + 8e-^ /2 , 

which, for r = 10 gives 1.0000012056138 . . . which is in error only in the 14th digit. Going 
further, we can write 



a 



l + 8e-" r/2 + 32e~ 7rr + 0(e- 37rr/2 ), (3) 



which, for r = 10 gives 19 significant digits. 

The preimage of the center of the rectangle will be on the imaginary axis, so write it 
as di. To find d in terms of a, consider rotating the rectangle by 180 degrees. This will 
correspond (under/) to a conformal automorphism of the upper half plane, i.e., a Mobius 
transformation 4>(z). Since the rotation interchanges opposite corners, <j) will interchange 
1 and —a and it will interchange 1 and a. The rotation interchanges the midpoints of the 
horizontal edges of the rectangle, so <fi interchanges and oo. So (f>(z) = a/z. The point 
di should be fixed by (ft, giving d 2 = a. 

For the random walk and the SAW in the half plane starting at id, it follows from 
([2]) and (fi]) that the (unnormalized) hitting density along the real axis is (x 2 + d 2 )~ h = 
(x 2 + a) . So letting pn denote the hitting density for a rectangle with the walk starting 
at the center, we have by ((2]) and Q, 

(x 2 + a)- b x\f'(x)\ b p R (f(z)). (4) 

We require the ratio of the integral of p(z) along a vertical edge to the integral along 
a horizontal edge, 

fo pn(.g/ 2 + iy)dy 

/-a/2 PR( x ) dx 

By a change of variable, setting u = f~ 1 (x) in the denominator and u = f~ 1 (a/2 + iy) in 
the numerator, this becomes 

Ji Pn(f(u))\f(u)\du 
S\p R (f(u))\f>(u)\du 

Recall that f'(u) = (1 — u 2 )~ l l 2 {a 2 — -u 2 ) -1 / 2 , so that the ratio of probabilities of a first 
hit on the vertical side to a first hit on the horizontal side, R(a, b) is 

, , r (n 2 + a)-\u 2 - l)( fe " 1 )/ 2 (a 2 - u^-^du 

R{a,b) = . (5) 

j\(u 2 + a)~ b (l - u 2 )i b - l )/ 2 {a 2 - u 2 )( b - l )/ 2 du 

4 Random walks. 

We first consider the random walk case, corresponding to b = 1. The integrals in ^ 
greatly simplify, giving 



^ arctan(- v /a) — arctan(l/y / a) 
2arctan(l/y / a) 
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It is straightforward to calculate the asymptotic expansion of both the numerator and 
denominator, and hence their ratio. In this way we find 

R(a, 1) = ^e~ wr/2 + O (e-" r ) . 

For r = 10 this evaluates to 3.83758797925134 . . . x 10 , which is correct to 6 significant 
digits. 

Proceeding to the second-order term, we obtain 

R(a, 1) = -e~^ 2 + + O (e'^l 2 ) . 

For r = 10 this evaluates to 3.8375894519594 . . . x 10~ 7 , which is correct to 13 significant 
digits. It is straightforward to obtain further terms should one so wish. 



5 The general case, including SAW. 

Unfortunately for 5 / 1 we cannot evaluate the integrals exactly, but as for the random 
walk case, asymptotics provides us with sufficient accuracy. Referring to equation ([5]), and 
noting that a is very close to 1 when the aspect ratio is 10, the denominator integral can 
be accurately approximated by 



\u 2 + l)-\l-u^du = ^M- 

1 2r U + 2) 



The numerator integral has a very narrow range of integration, so we can approximate 
the factor u 2 + a by 2, the factor (u 2 — 1) by 2(u — 1) and the factor a 2 — u 2 by 2(a — u). 
Thus the numerator can be approximated by 



l)( 6 -D/2 (a _ u )(6-l)/2 du . 



If we set u = 1 + t(a — 1), this becomes 

ha - If /V " t)\ {b - 1)/2 dt = ^'fl^ ^ ~ l ) h - 

1 JO T \ 1 +2) 

Combining the above results, and using the result derived above, a — 1 ~ 8e _7rr//2 , we 
find, asymptotically 

R(r, b) « H^ L e~ nbr/2 . (6) 

r(i + |)r(|) 



For SAW, 6 = 5/8, and @ reduces to 

R(r, 5/8) w 1.2263431442e- 57rr/16 . 
For r = 10 this gives ^(10,5/8) ps 6.6824528 x 10" 5 . 
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In order to get an estimate of the accuracy of this result, we calculated the relevant 
integrals numerically, which is quite straightforward in either Maple or Mathematica. In 
this way we found the more accurate estimate R(10, 5/8) ~ 6.682989935 x 1CP 5 . 

A rather elaborate calculation gives the next term in the asymptotic expansion of the 
aspect ratio as 



R(r,b) 



i26+l 



A 



-bnr/2 



1 + 



A2 b+1 

2 



-bwr/2 



+ 4(6-l + 2A)e-^ r/2 + 0(e 



for < b < 1, where A = ( I . Evaluating this for r = 10, b = §, we find 



£(10, 5/8) = 0.00006682989679 • • • 

which is in error by 2 in the 8 th significant digit, whereas the leading order term gives an 
estimate which is in error by 6 in the 5 th significant digit 

We see that the corresponding ratio for the scaling limit of SAWs in a 10 x 1 rectangle 
is some two orders of magnitude larger than the corresponding result for Brownian motion. 
This higher probability is intuitively what one might expect, as SAWs tend to move in a 
straight trajectory rather more than do random walks. 
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